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ABSTRACT 


As  part  of  the  Aircraft  Battle  Damage  Repair 
(ABDR)  Project,  CNA  analysts  performed  research  into 
finding  an  appropriate  model  for  determining  ABDR 
spares  requirements.  The  analysis  focuses  on  the  rela¬ 
tive  performance  of  four  spares  requirements  models, 
given  uncertainties  associated  with  predicting  battle 
damage  rates. 
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SUMMARY  AND  RESULTS  OF  ANALYSIS 


INTRODUCTION 

%  As  part  of  the  Aircraft  Battle  Damage  Repair  (ABDR)  Study,  CNA  analysts 

performed  research  into  finding  an  appropriate  model  for  determining  ABDR 
material  requirements.  This  research  memorandum  documents  the  analysis  and 
recommends  a  specific  model  for  determining  ABDR  spues  requirements. 

This  section  of  the  paper  contains  a  nontechnical  discussion  of  the  analysis, 
a  summary  of  the  results,  and  the  conclusion  of  the  analysis.  The  next  section 
contains  technical  details,  including  a  discussion  of  the  mathematics  behind  the 
analysis  and  precise  definitions  of  the  spares  models  and  probability  distributions 
used  in  the  analysis.  Detailed  results  of  the  analysis  are  presented  in  appendix  A. 
Two  of  the  spares  models  considered  in  this  analysis  had  a  surprising  relationship; 
Appendix  B  contains  a  detailed  example  that  helps  to  illustrate  the  reasons  for 
this.  A  listing  of  the  major  computer  program  used  in  the  analysis  is  contained 
in  appendix  C. 


BACKGROUND 

The  overall  modeling  process  for  ABDR  material  requirements  may  be  broken 
into  two  parts:  (1)  the  development  data  bases  for  predicting  battle  damage 
rates,  and  (2)  development  of  a  spares  model  to  develop  ABDR  spares  lists  from 
these  data  bases.  Current  Navy  planning  calls  for  these  data  bases  for  predicting 
battle  damage  rates  to  be  developed  by  the  Naval  Weapons  Center  (NWC)  and 
the  Naval  Air  Development  Center  (NADC).  Essentially,  NWC  vulnerability  and 
susceptibility  models  will  be  used  to  model  aircraft  damage;  the  ABDR  material 
required  to  repair  each  type  of  damage  incident  is  determined  in  conjunction  with 
NADC  data  bases.  The  end  result  of  the  process  will  probably  be  a  collection 
of  component  or  part-level  battle  damage  rates.  Following  the  development  of  a 
data  base  for  predicting  battle  damage  rates,  NADC  will  use  a  spares  model  to 
develop  ABDR  allowance  lists. 

CNA’s  research  effort  focused  only  on  the  spares  model  and  did  not  consider 
the  methods  used  to  develop  the  data  base  for  predicting  battle  damage  rates. 
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The  objective  of  the  research  effort  was  to  find  a  spares  model  appropriate  for 
determining  ABDR  material  requirements. 

There  are  at  least  two  striking  differences  between  the  types  of  failures  expe¬ 
rienced  normally  in  peacetime  and  those  associated  with  battle  damage:  uncer¬ 
tainty  and  dependence.  4 

The  peacetime  failure  rates  from  which  Aviation  Consolidated  Allowance 
Lists  (AVCALs)  are  developed  are  based  on  actual  usage  data.  In  determining 
ABDR  requirements,  planners  must  use  predicted  battle  damage  rates  rather 
than  peacetime  usage  data;  of  course,  there  are  not  any  relevant,  actual  battle 
damage  data  for  use  in  this  process.  The  Navy’s  plans  to  use  vulnerability  and 
susceptibility  models  to  develop  predicted  battle  damage  rates  seems  reasonable; 
indeed,  an  alternative  method  for  developing  predicted  battle  damage  rates  is 
not  evident.  However,  regardless  of  how  detailed  and  accurate  the  models  are 
from  an  engineering  standpoint,  the  damage-rate  predictions  will  be  based  in 
large  part  on  the  following  information: 

•  Educated  guesses  for  model  parameters,  particularly  descriptions  of  threats, 
tactics,  and  operational  objectives 

•  Intelligence  information 

•  Historical  data  from  conflicts  in  which  weapons,  tactics,  and  threats  were 
different  from  those  likely  to  be  experienced  in  a  future  conflict. 


Hence,  there  is  more  uncertainty  associated  with  battle  damage  rates  than  with 
empirical,  peacetime  usage  data  used  to  develop  AVCAL  requirements.  This 
uncertainty  should  be  considered  when  selecting  a  model  for  ABDR  material 
requirements. 

Peacetime  spares  models  assume  independence  of  failures.  That  is,  given  two 
distinct  items  on  the  aircraft,  the  probability  that  item  A  fails  during  a  flight  is 
independent  of  the  probability  that  item  B  fails  during  the  flight.  In  general,  the 
assumption  of  the  independence  of  battle  damage  “failures”  is  not  easy  to  accept. 
Indeed,  it  is  clear  that  the  physical  proximity  of  two  items  is  important,  and  if 
two  items  are  located  close  together,  they  will  tend  to  be  damaged  together. 

Consider  the  following  example  of  the  effect  of  the  independence  assumption 
in  the  battle  damage  situation.  Assume  two  parts,  A  and  B,  are  both  required  to 
have  a  flyable  aircraft.  Also  assume  that  parts  A  and  B  have  identical  predicted 
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battle  damage  rates,  and  part  A  costs  ten  times  as  much  as  part  B.  A  readiness- 
based  spares  model,  analogous  to  the  readiness-based  AVCAL  model  currently 
planned  for  use  in  developing  F-18  AVCALs,  trades  off  cost  against  improvements 
in  readiness  when  computing  allowance  levels.  This  type  of  spares  model  also 
assumes  the  independence  of  failures.  The  readiness-based  model  would  tend  to 
stock  part  B  to  a  higher  level  than  part  A.  Just  for  the  sake  of  this  example, 
suppose  ten  items  of  part  B  are  stocked  and  four  items  of  part  A  are  stocked. 
Assume  that  parts  A  and  B  are  located  next  to  each  other  in  the  aircraft,  so  that 
part  A  and  part  B  are  always  damaged  together.  In  this  case,  six  extra  items  of 
part  B  are  useless,  because  any  aircraft  needing  the  replacement  of  part  B  also 
requires  the  replacement  of  part  A. 

Clearly,  the  assumption  of  independence  in  a  peacetime  spares  model  must 
be  carefully  considered  if  the  model  is  adapted  for  ABDR  spares. 

Any  model  for  developing  spare  parts  requirements  must  use  predicted  battle 
damage  rates.  As  discussed  above,  actual  battle  damage  rates  may  be  different 
from  the  predicted  rates;  moreover,  the  actual  rates  may  not  satisfy  assumptions 
of  independence  built  into  certain  spares  models.  Most  spares  models  will  do 
well  in  determining  requirements  when  the  predicted  and  actual  damage  rates 
are  equal,  but  what  happens  if  the  predicted  rates  and  actual  damages  rates  are 
different?  The  following  example,  depicted  in  figure  1,  illustrates  the  idea  that 
ABDR  sparing  models  should  have  the  flexibility  to  provide  support  under  the 
uncertainty  assumption  about  the  battle  damage  rate  planning  factors.  Suppose 
two  different  models  are  used  to  form  ABDR  spare  parts  inventories,  A  and  B 
respectively,  from  the  same  set  of  predicted  damage  rates.  Suppose  stockpile 
A  performs  much  better  than  stockpile  B  when  the  actual  damage  rates  and 
the  predicted  damage  rates  are  equal,  but  as  actual  damage  rates  and  predicted 
damage  rates  become  further  apart,  stockpile  B’s  performance  becomes  better 
than  stockpile  A’s.  Stockpile  A,  optimized  for  the  situation  in  which  predicted 
and  actual  damage  rates  are  equal,  is  not  flexible  enough  to  cope  well  with 
different  damage  rates.  Stockpile  A  performs  well  only  in  a  specific,  narrow 
setting.  Stockpile  B,  designed  for  flexibility,  performs  reasonably  well  in  a  wide 
range  of  situations  but  is  not  the  optimal  model  when  predicted  and  actual 
damage  rates  are  equal. 
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Effectiveness  of  stockpile 


* 


RG.1:  FLEXIBILITY  VERSUS  OPTIMALITY 
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NONTECHNICAL  DEFINITIONS  OF  SPARES  MODELS 


Four  alternative  spares  models  were  considered  in  this  analysis:  the  DBS- 
Part,  RBS-Part,  DBS-Damage  Incident,  and  RBS-Damage  Incident  models.1 
Given  a  list  of  candidate  parts  for  inclusion  in  a  spare  parts  stockpile,  each 
spares  model  is  designed  to  develop  a  stock  level  for  each  candidate  part  suffi¬ 
cient  to  cover  demands  over  a  user-specified  period  of  time  under  the  assumption 
of  no  resupply  or  repair  of  the  damaged  part. 

The  list  of  candidate  parts  does  not  necessarily  include  all  parts  on  the  air¬ 
craft.  For  example,  the  Navy  may  decide  to  to  include  only  items  that  are  not 
normally  required  to  support  peacetime  flight  operations;  in  particular,  AVCAL 
items  would  probably  not  be  included  in  an  ABDR  spares  stockpile.  In  the 
following  analysis,  “damage”  refers  to  the  need  to  replace  a  candidate  part  on 
the  aircraft.  Battle  damage  not  leading  to  a  demand  for  a  replacement  part  is 
treated  as  “no  damage”  by  these  models.  In  particular,  damage  requiring  the 
replacement  of  “noncandidate”  parts  is  treated  as  “no  damage”  by  the  spares 
models. 

The  DBS-Part  model  is  a  direct  adaptation  of  the  demand-based  Aviation 
Supply  Office  (ASO)  spares  model  currently  used  to  develop  AVCALs  for  most 
aircraft  types.  The  stock  level  of  a  part  under  this  rule  depends  only  on  the 
predicted  damage  rate  for  the  part.  Assuming  damage  incidents  are  described 
by  a  binomial  distribution,  just  enough  parts  are  stocked  to  ensure  that  the 
probability  of  a  stockout  is  no  more  than  10  percent. 

The  RBS-Part  model  is  analogous  to  the  readiness-based  spares  model  used 
by  ASO  to  develop  AVCAL  lists  for  the  F-18.  The  stock  level  of  a  part  depends 
on  the  unit  cost  of  the  part  as  well  as  the  predicted  damage  rate  of  the  part. 
Essentially,  the  part  stockpile  is  built  part  by  part,  where  at  each  step,  one  unit 
of  the  item  offering  the  largest  gain  in  “readiness”  per  unit  cost  is  added.  The 
process  is  stopped  when  a  specified  cost  goal  is  reached.  In  this  analysis,  when 
compared  to  the  DBS  models,  the  total  cost  for  the  stockpile  produced  by  RBS- 
Part  was  constrained  to  be  just  below  the  cost  of  the  corresponding  DBS-Part 
stockpile.  In  other  parts  of  the  analysis,  other  cost  goals  were  used;  in  each  case, 
the  cost  goal  used  is  clearly  indicated  in  this  report. 

1.  DBS  stands  for  demand-based  sparing,  and  RBS  stands  for  readiness- based  sparing.  DBS-P, 
RBS-P,  DBS-DI,  and  RBS-DI  will  be  used  to  denote  the  DBS  part,  RBS  part,  DBS  damage 
incident,  and  RBS  damage  incident  models,  respectively. 
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Both  of  the  spares  models  described  above  assume  independence  between  part 
damage  rates.  As  discussed  earlier,  this  assumption  does  not  seem  reasonable 
when  applied  to  battle  damage.  The  “damage  incident”  models  are  designed  to 
replace  the  unacceptable  assumption  of  independence  between  part  damage  with 
the  more  acceptable  assumption  of  independence  between  damage  incidents. 

The  DBS-Damage  Incident  model  is  a  direct  adaptation  of  the  DBS-Part 
model,  where  damage  incidents  play  the  role  of  “parts”  in  the  DBS-Part  model. 
In  effect,  “bundles”  of  parts,  as  defined  by  the  dependence  among  parts  in  battle 
damage  incidents,  would  be  carried  in  the  ABDR  part  stockpile.  For  example, 
suppose  there  are  two  damage  incidents:  incident  A  requires  parts  1  and  2,  while 
incident  B  requires  parts  2  and  3  for  repairs.  Suppose  the  probability  of  expe¬ 
riencing  more  than  four  of  incident  A  is  less  than  10  percent,  and  experiencing 
more  than  three  of  incident  B  is  less  than  10  percent.  Then  the  “damage  incident 
stockpile”  is  four  of  incident  A  and  three  of  incident  B.  In  terms  of  parts,  this 
corresponds  to  the  following  stockpile:  four  of  part  1 ,  seven  of  part  2,  and  three 
of  part  3. 

In  an  analogous  manner,  the  RBS-Damage  Incident  model  is  adapted  from 
the  RBS-Part  model. 


COMPARISON  OF  THE  ALTERNATIVE  SPARES 
MODELS 

This  analysis  is  based  on  the  premise  that  predicted  damage  rates  and  actual 
damage  rates  are  likely  to  be  different,  and  a  spares  model  offering  flexibility  (e.g., 
model  B  in  figure  1)  is  preferable  to  a  spares  model  offering  narrow  optimality 
(e.g.,  model  A  in  figure  1).  The  flexibility  of  the  four  alternative  spares  models 
was  investigated  over  a  range  of  three  possible  relationships  between  predicted 
and  actual  damage  rates: 


•  Full  confidence:  Actual  damage  rates  are  the  same  as  the  predicted  damage 
rates. 

•  No  confidence:  No  relationship  or  correlation  exists  between  predicted  and 
actual  damage  rates. 

•  Some  confidence:  Predicted  rates  are  close  to  the  actual  damage  rates. 
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A  precise  definition  of  the  “confidence  factor”  parameter  used  in  this  analysis  is 
given  in  the  technical  section  of  this  paper. 

,  The  models  were  compared  in  terms  of  a  random  variable  CANN,  roughly  cor¬ 

responding  to  the  number  of  cannibalized  aircraft  needed  to  meet  the  deficiencies 
in  stockpiles  of  ABDR  spares.  A  simple  illustration  of  the  way  CANN  is  calcu- 
t  lated  is  given  in  table  1. 

TABLE  1 

EXAMPLES  OF  CALCULATING  CANN 
_ Part  1 _ _  _ Part  2 _ 

Example 

number  Stocked  Demanded  Stocked  Demanded  CANN 

1  3  5  7  6  2 

2  3  2  7  5  0 

3  3  7  7  12  5 

Initial  work  was  based  on  a  list  of  23  major  structural  assemblies  for  the  F-14 
that  the  Naval  Air  Systems  Command  (NAVAIR)  is  considering  as  candidates 
for  ABDR  spares  stockpiles.  A  set  of  30  damage  incidents  was  hypothesized, 
and  the  performance  of  the  four  spares  models  in  this  situation  was  analyzed. 
At  later  stages  of  the  analysis,  more  general  situations  were  considered  in  which 
both  the  cost  of  ABDR  spares  and  the  list  of  parts  associated  with  each  damage 
incident  were  described  by  probability  distribution  functions.  The  generalized 
cost  distribution  used  was  derived  from  cost  data  for  major  aircraft  assemblies 
provided  by  NAVAIR.  The  generalized  distributions  of  part-to-damage  incidents 
were  designed  to  allow  the  parameterization  of  the  level  of  dependence  between 
part  damage  rates. 

The  numbers  of  puts  (23)  and  damage  incidents  (30)  used  in  the  initial  work 

♦  on  F-14  major  assemblies  were  retained  throughout  the  analysis.  Because  of 
computer  limitations,  it  was  not  feasible  to  extend  the  analysis  to  substantially 
larger  numbers  of  parts  and  damage  incidents.  This  was  not  judged  to  severely 

*  limit  the  generality  of  the  analysis,  because  a  wide  range  of  damage  rates  and 
distributions  of  part-to-damage  incidents  were  considered  in  the  comparison  of 
the  spares  models. 
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A  sample  of  results  is  contained  in  figure  2,  and  tables  2  and  3.1  Figure  2 
compares  the  probability  distributions  of  CANN.  Tables  2  and  3  show  a  more 
direct  comparison  between  the  spares  models.  Model  A  is  said  to  be  better  than 
model  B  if  model  A  leads  to  fewer  Cannibalized  aircraft”  than  model  B;  that 
is,  if  CANNfl  —  CANN,*  >  0.  These  sample  results  correspond  to  a  stockpile 
of  major  aircraft  assemblies  and  are  based  on  a  support  period  of  3,000  sorties 
without  resupply,  an  overall  damage  rate  of  approximately  0.008  per  sortie,  and 
“some  confidence”  in  the  predicted  damage  rates.  In  figure  2  and  table  2,  the  cost 
of  the  RBS  stockpile  is  limited  to  the  cost  of  the  corresponding  DBS  stockpile; 
in  this  comparison,  the  cost  of  a  Damage  Incident  stockpile  is  higher  than  the 
cost  of  the  Part  stockpile.  Table  3  shows  an  equal-cost  comparison  between  the 
RBS  models. 


TABLE  2 


SAMPLE  RESULTS:  COMPARISON  OF 
FOUR  SPARES  MODELS 


The  “probability  Model  A  performs  better  than 
Model  B”  is  Prob{CANNB  -  CANN*  >  0}. 

Parameters  used  in  this  comparison 
RBS-P  cost  goal:  DBS-P  cost 
RBS-DI  cost  goal:  DBS-DI  cost 
Some  confidence  in  predicted  damage  rates 

DBS-P  versus  RBS-P  DBS-P  versus  DBS-DI 


Prob.  DBS-P  is  better:  0.11  Prob.  DBS-P  is  better:  0.01 
Prob.  RBS-P  is  better:  0.41  Prob.  DBS-DI  is  better:  0.29 


RBS-P  versus  DBS-DI 

Prob.  DBS-P  is  better:  0.06 
Prob.  RBS-DI  is  better:  0.48 


RBS-P  versus  DBS-DI 


Prob.  RBS-P  is  better:  0.28 
Prob.  DBS-DI  is  better:  0.21 


DBS-P  versus  RBS-DI  RBS-DI  versus  RBS-DI 

Prob.  RBS-P  is  better:  0.06  Prob.  DBS-DI  is  better:  0.12 
Prob.  RBS-DI  is  better:  0.24  Prob.  RBS-DI  is  better:  0.34 

1.  Tables  2  and  3  are  taken  bom  “series  5"  of  appendix  A. 
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(CANN  >  N) 


* 


N 


FIG.  2:  SAMPLE  RESULTS:  PROBABILITY-DISTRIBUTION 
FUNCTIONS  OF  CANN 
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TABLE  3 


SAMPLE  RESULTS:  COMPARISON  OF  * 

EQUAL-COST  GOALS  OF  RBS  MODELS 

The  “probability  Model  A  performs  better  than  Model  B”  J 

is  Prob{CANNB  -  CANN^  >  0}. 

Parameters  used  in  this  comparison 
RBS-P  cost  goal:  slightly  higher  than  DBS-DI  cost 
RBS-DI  cost  goal:  slightly  higher  than  DBS-DI  cost 
Some  confidence  in  predicted  damage  rates 

RBS-P  versus  RBS-DI 


Prob.  RBS-P  is  better:  0.23 
Prob.  RBS-DI  is  better:  0.05 


The  relationship  between  the  RBS  and  DBS  models,  illustrated  in  figure  2 
and  table  2,  held  consistently  among  the  situations  studied.  The  RBS  (Part  or 
Damage  Incident)  approach  was  consistently  better  over  all  three  levels  of  con¬ 
fidence  in  predicted  battle  damage  rates  than  the  corresponding  DBS  approach. 
The  results  indicate  that  the  readiness-based  approach  produces  a  better  spares 
stockpile  than  the  demand-based  approach. 

Recall  that  in  the  comparison  in  table  2,  the  RBS-Part  cost  goal  was  set 
so  that  the  cost  of  the  RBS-Part  stockpile  was  slightly  below  the  cost  of  the 
DBS-Part  stockpile.  It  is  always  true  that  the  cost  of  the  DBS-Damage  Incident 
stockpile  will  be  greater  than  the  cost  of  the  DBS-Part  stockpile.  Hence,  in  the 
table  2  comparison,  the  DBS-Damage  Incident  stockpile  is  higher  in  cost  than 
the  RBS-Part  stockpile.  The  relationship  between  the  DBS-Damage  Incident  and 
RBS-Part  models  illustrated  in  table  2  is  surprising  in  view  of  this  relative  cost 
relationship.  The  RBS-Part  model  generally  performed  better  than  the  DBS- 
Damage  Incident  model  when  predicted  damage  rates  were  not  good  predictions 
of  actual  damage  rates;  however,  the  opposite  was  true  when  the  predicted  dam¬ 
age  rates  were  good  predictions  of  the  actual  damage  rates.  Hence,  the  RBS-Part 
stockpile  generally  performed  better  than  the  more  expensive  DBS-Damage  In¬ 
cident  stockpile.  This  result  emphasizes  the  flexibility  of  the  readiness-based 
approach. 
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Recall  that  in  table  2,  the  RBS-Part  and  RBS-Damage  Incident  models  had 
different  cost  goals,  making  it  difficult  to  compare  these  two  models.  To  differ¬ 
entiate  their  relative  value,  analysts  made  a  number  of  equal-cost  comparisons 
between  the  RBS-Part  and  RBS-Damage  Incident  model.  Because  of  the  high 
cost  of  some  parts,  even  when  the  two  models  are  run  to  equal-cost  goals,  there 
will  probably  be  some  (relatively  small)  difference  in  the  cost  of  the  stockpiles. 

In  the  equal-cost  comparison  illustrated  in  table  3,  the  cost  goals  were  set  so 
that  so  that  the  RBS-Part  stockpile  cost  would,  on  average,  fall  slightly  below  the 
cost  of  the  RBS-Damage  Incident  model.  In  effect,  this  gives  a  slight  advantage  to 
the  RBS-Damage  Incident  model.  Table  3  contains  typical  results;  the  RBS-Part 
model  generally  performed  better  than  the  RBS-Damage  Incident  model  in  equal- 
cost  comparisons.  In  fact,  the  RBS-Damage  Incident  had  better  performance 
only  in  a  few  cases:  a  subset  of  the  comparisons  in  which  the  damage  rates  were 
high  and  there  was  little  correlation  between  predicted  and  actual  damage  rates. 

The  excellent  performance  of  the  RBS-Part  model  was  unexpected;  it  was 
anticipated  that  the  independence  assumption  in  the  RBS-Part  model  would  re¬ 
sult  in  poor  performance  relative  to  the  RBS-Damage  Incident  model.  From  the 
way  the  cost  goals  of  table  3  were  set,  the  RBS-Part  stockpile  costs  slightly  less 
on  average;  hence,  the  better  performance  is  not  due  to  a  difference  in  cost.  Rel¬ 
ative  to  a  DBS  model,  the  RBS  model  will  reduce  the  stock  level  of  expensive 
items  and  use  the  savings  to  increase  stock  levels  of  less  expensive  parts.  In  this 
way,  the  probability  of  experiencing  stockouts  of  inexpensive  items  is  “greatly” 
reduced  at  the  expense  of  “slightly  increasing”  the  probability  of  stockouts  for 
the  expensive  items.  The  RBS-Damage  Incident  model  appears  to  carry  this  ten¬ 
dency  of  trading  off  expensive  items  for  inexpensive  items  a  bit  too  far.  Relative 
to  the  RBS-Part  model,  the  RBS-Damage  Incident  model  appears  to  experience 
slightly  more  stockouts  for  expensive  items  and  slightly  fewer  stockouts  for  inex¬ 
pensive  items.  However,  the  RBS-Part  model  also  does  well  in  meeting  demands 
for  the  inexpensive  items.  The  overall  effect  is  that  the  RBS-Part  model  gen¬ 
erally  performs  slightly  better  than  the  RBS-Damage  Incident  model.  A  more 
detailed  example  is  provided  in  appendix  B. 

Appendix  A  contains  ProbfCANN*  -  CANNb  >  0},  as  in  tables  2  and  3,  for 
the  different  situations  considered  in  the  analysis.  While  graphs  corresponding 
to  figure  2  were  also  considered  in  the  analysis,  they  did  not  add  information  and 
are  not  included  in  appendix  A. 
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CONCLUSION 


Among  the  four  spares  models  considered  in  this  analysis,  the  RBS  models 
provide  better  performance  than  the  corresponding  DBS  models.  When  RBS- 
Part  and  RBS-Damage  Incident  are  used  to  equal-cost  goals,  the  RBS-Part  model 
performs  better  overall.  In  addition,  RBS-Part  is  an  easier  model  to  use  because 
damage  rates  and  similar  data  are  processed  at  the  part  level  rather  than  the 
damage-incident  level.  The  analysis  indicates  that  the  RBS-Part  model  should 
be  adopted  for  determining  ABDR  material  requirements. 
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TECHNICAL  DETAILS 


In  this  section,  a  number  of  technical  details  are  given,  including  a  discussion 
of  the  mathematics  used  in  the  analysis  and  precise  definitions  of  the  spares 
models  and  distribution  functions  used  in  the  analysis. 

MATHEMATICS 

The  underlying  events  in  the  mathematical  formulation  of  the  ABDR  spares 
model  problem  are  damage  incidents.  In  terms  of  the  NWC  models  used  to  study 
aircraft  vulnerability  and  susceptibility,  each  damage  incident  may  correspond  to 
a  specific  “shotline”  in  a  vulnerable  area  model,  or  to  a  specific  missile  endgame 
model  result.  (See  [1]  for  a  discussion  of  these  models.)  In  real-world  terms,  the 
set  of  damage  incidents  is  the  set  of  possible  outcomes  of  a  combat  sortie.  For 
mathematical  convenience,  the  outcome  in  which  no  parts  are  required,  either 
because  no  damage  is  experienced  or  none  of  the  candidate  parts  are  needed  for 
repairs,  will  always  be  included  in  the  set  of  damage  incidents.  One  assumption 
made  with  respect  to  damage  incidents  is  that  the  damage  incident  experienced 
on  each  sortie  is  independent  of  the  damage  incidents  experienced  on  other  sor¬ 
ties. 

Assume  there  is  a  set  of  np  parts  that  are  candidates  for  an  ABDR  stockpile, 
and  assume  that  there  is  a  set  of  nj  distinct  damage  incidents.  The  following 
notation  is  used  throughout  this  section  to  denote  ABDR-related  factors  over 
some  period  of  time  T  encompassing  N  combat  sorties: 


•  The  vectors  0  and  1  denote  (0, 0,  •  •  • ,  0)  and  (1, 1,  •  ■  • ,  1),  respectively.  Vec¬ 
tor  addition,  scalar  multiplication,  and  “dot  product”  (denoted  •)  have  the 
usual  meaning.  Two  vectors  satisfy  v  <  w  if  the  inequality  (tf)<  <  (u>), 
holds  for  the  i  components  of  the  vectors. 

•  The  vector  rp  describes  predicted  damage  rates.  The  i-th  component  de¬ 
notes  the  probability  that  damage  incident  i  will  be  experienced  on  a  sortie. 
Recall  that  the  no-damage  incident  is  included  with  appropriate  rp  and  note 
that  rp  •  1  =  1. 
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•  The  vector  r„  describes  actual  damage  rates.  These  are  damage  rates  that 
will  be  experienced  in  combat  and  may  be  different  from  the  predicted 
battle  damage  rates. 

•  The  vector  k  describes  the  number  of  each  type  of  damage  incident  experi- 
enced  over  a  period  of  time  T.  The  i-th  component  (k)i  denotes  the  number 
of  times  damage  incident  i  was  experienced.  The  probability  distribution 
of  k  after  N  combat  sorties  is  the  multinomial  distribution: 

Prob{fc  =  x  after  N  sorties}  =  M(x,r^,JV) 


•  The  vector  D{  specifies  the  list  of  parts  required  to  repair  damage  incident 

i.  For  example,  if  D{  =  (1, 0, 1,  •  •  •)»  then  one  of  part  number  1,  one  of  part 

number  3,  and  so  on,  are  needed  to  repair  damage  incident  i.  This  analysis 
— *  -♦ 

assumes  that  Di  <  1 

•  The  vector  d  describes  the  number  of  each  type  of  part  demanded  over 

the  time  period  T.  Given  the  vector  k  describing  the  damage  incidents 
experienced  this  period,  (d)i  =  Note  that  d  and  are 

dependent  random  vectors;  and  d  also  depends  on  the  number  of  combat 
sorties  N,  and  the  part-to-damage  incident  vectors  D. 

•  The  vector  s  describes  stock  levels.  The  i-th  component  (s)i  denotes  the 
stock  level  of  item  i  in  the  ABDR  stockpile.  Similarly,  c  describes  the 
costs  of  the  ABDR  parts.  The  stock  level  produced  by  a  spares  model  is 
a  function  of  the  predicted  damage  rates,  the  number  of  combat  sorties 
N  to  be  supported  by  the  stockpile,  cost,  and  the  part-to-damage-incident 
vectors  D ;  in  functional  notation,  a  =  a{r^,N,c,D). 


Measure  of  Effectiveness  Used  To  Compare  Spares  Models 

The  models  for  spare  parts  requirements  were  compared  in  terms  of  a  random 
variable  CANN,  computed  by  comparing  the  list  of  parts  demanded  to  the  list 
of  parts  stocked.  Let  (s)<  be  the  stock  level  of  part  i,  and  let  (d)<  be  the  demand 
for  part  i.  Assuming  no  resupply,  max{0,  (d),  -  (s};}  is  the  number  of  stockouts 
for  part  i.  The  random  variable 


CANN  =  max{max{0,  (d),  —  (s%}|  all  i} 
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is  analogous  to  the  number  of  cannibalized  aircraft  required  to  fill  the  stockouts. 
The  notations  CANN.*  and  CANN b  are  used  to  represent  CANN  under  the 
sparing  models  A  and  B,  respectively.  The  random  variable  CANN.*  —  CANN# 
is  used  to  directly  compare  spares  model  A  with  spares  model  B. 

It  is  convenient  to  think  of  CANN  as  a  function  of  the  stock  level,  demand  for 
parts,  predicted  and  actual  damage  rates,  the  number  of  combat  sorties  N,  cost, 
and  the  part-to-damage-incident  vectors  D;  in  functional  notation,  CANN  = 
CANN(i*,d,ra,rp,  iV,c,.D).  Also,  it  is  convenient  to  think  of  the  variables  of 
interest  as  belonging  to  a  probability  space  X  =  {(i*,  d,  r^,  r£,  N,  c,  D)}  equipped 
with  a  probability  measure  fi.  In  this  space,  N,  c,  and  D  are  all  independent 
of  other  variables;  s  is  completely  determined  by  7*p,  c,  and  JV;  and  there  are 
dependencies  between  d,  r^,  and  r^.  This  allows  certain  quantities  of  interest  in 
the  analysis,  in  particular  Prob{CANN*  —  CANNb  >  0},  to  be  represented  as 
an  integral 


for  some  specific  function  /.  If  /  is  defined  by 


f{s,d,ra,rp,N,c,D) 


1  if  CANN(-)x  -  CANNb(-)  >  0 

0  otherwise 


then  I  =  Prob{CANN,*  -  CANNb  >  0}. 

The  integral  I  introduced  above  is  difficult  to  calculate  because  the  number 
of  numerical  calculations  needed  to  compute  the  exact  value  of  the  integral  grows 
rapidly  with  the  number  of  damage  incidents  and  parts  considered.  Stochastic 
integration  methods  were  selected  as  the  easiest  way  of  estimating  the  integral. 
Justification  of  the  techniques  used  are  given  below.  A  good  general  reference 
on  stochastic  integration  is  [2]. 

Estimates  of  the  quantities  Prob{CANN*  —  CANNb  >  0}  were  made  over 
a  variety  of  damage  rates,  costs,  and  part-to-damage-incident  tables,  each  com¬ 
bination  giving  rise  to  a  distinct  probability  measure  fi.  The  parameters  and 
distribution  functions  used  to  describe  damage  rates,  costs,  and  part-to-damage- 
incident  tables  are  given  later  in  this  section.  The  estimates  are  reported  in 
appendix  A. 
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Stochastic  Integration 


The  integrals  considered  in  this  analysis  are  of  the  form  I  =  fx  fdfi ,  where 
/  :  X  — *  {0,1}  and  /i  is  a  probability  measure  on  X.  The  statistic  0„  is  defined 
by  0n  =  1  f(xi ),  where  {ajj  :  t  =  1,  •  •  • ,  n}  is  an  independent  random  sample 

of  size  n  from  X. 

The  expected  value  of  the  statistic  dn  is 

E(K)  =  fxfd^  =  I  , 

and  the  variance  of  the  statistic  is 

V{K)  =  -([  fi»-P)  . 

n  J  x 

By  Chebyshev’s  inequality,  and  using  the  facts  that  I  —  I2  <  1/4  on  [0,1]  and 

/2  =  /, 

Prob{|*n-/|>e}<iF(0n) 

=  hI~I')^ik  ■ 

This  allows  one  to  determine  the  number  of  samples  required  to  give  a  specified 
level  of  accuracy  to  the  estimate  0n.  The  sample  size  used  in  all  estimates  reported 
in  this  analysis  was  n  =  500,  ensuring  for  example  that 

Prob{|0n  -  I\  >  0.1}  <  0.05; 

that  is,  95  percent  confidence  that  the  estimate  is  within  0.1  of  the  true  value. 
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Numerical  Calculations 


Two  separate  computer  programs  were  written  to  carry  out  the  stochastic 
integration  technique  described  above.  The  first  program  generates  500  random 
samples  from  the  probability  space  X.  For  each  of  the  random  samples,  four 
values  of  CANN,  one  for  each  of  the  four  spares  models,  are  calculated  and 
written  to  an  intermediate  file.  A  listing  of  this  program,  corresponding  to  the 
most  general  definition  of  the  probability  space  X  considered,  is  included  in 
appendix  C.  The  second  computer  program  reads  values  of  CANN  from  the 
intermediate  file,  then  computes  and  averages  the  corresponding  values  of  the 
function  /. 

An  outline  of  the  procedure  used  to  generate  a  single  random  sample  from 
the  probability  space  X  =  {(i*,d,ra,rp,  N,  c,  D)}  is  given  below: 


•  N:  In  all  of  the  situations  considered,  the  number  of  sorties  N  was  held 
constant. 

•  c:  In  some  situations  considered,  the  cost  vector  c  was  constant.  In  other 
situations,  a  random  vector  c  was  drawn  from  the  (independent)  probability 
distribution  for  cost. 

•  D:  In  some  situations  considered,  the  matrix  describing  the  list  of  parts 
needed  to  repair  a  damage  incident  was  constant.  In  other  situations,  a  ran¬ 
dom  matrix  D  was  drawn  from  the  (independent)  probability  distribution 
describing  part-to-damage-incident  relationships. 

•  f pi  A  random  vector  describing  the  predicted  damage  rates  was  drawn  from 
the  probability  distribution  for  r"p. 

•  r„:  The  predicted  and  actual  damage  rates  are  not  necessarily  independent. 
After  a  specific  predicted  damage  rate  vector  rp  was  drawn,  an  actual  dam¬ 
age  rate  vector  was  drawn  from  the  conditional  distribution  function  of 
ra  given  r^. 

•  d :  The  random  vector  k  describing  damage-incident  experience  is  depen¬ 
dent  on  the  actual  damage  rate  r„.  After  a  specific  actual  damage- rate 
vector  r"a  was  drawn,  a  damage-experience  vector  k  was  drawn  from  the 
conditional  distribution  function  of  k,  given  r„  and  N.  The  list  of  parts 
required,  d,  was  then  computed  directly  from  k  using  the  part-to-damage- 
incident  matrix  D. 
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•  i*:  The  vector  of  stock  levels  s  was  computed  directly  from  the  cost  vector 
c,  N,  and  the  vector  of  predicted  damage  rates  rp  using  the  appropriate 
spares  model  algorithm.  » 


DEFINITION  OF  SPARES  MODELS 

Four  different  spares  models  were  evaluated;  each  model  is  defined  below. 
In  the  following, 


B(k,p,n)  =  Y,  (j)  p’O-  ~P)n  3 
j= o  '  ' 

denotes  the  binomial  distribution  function.  Each  model  assumes  that  the  ABDR 
spares  are  to  cover  demands  over  a  fixed  period  of  time,  that  is,  a  fixed  number 
of  combat  sorties,  with  no  resupply  or  repair  of  damaged  parts. 

The  DBS-Part  model  is  a  direct  adaptation  of  the  supply-oriented  model 
traditionally  used  in  developing  AVCALs.  The  stock  level  of  part  i  is  the  first 
integer  k  satisfying  B(k,ri,N)  >  0.9,  where  r,  =  is  the  predicted 

damage  rate  for  part  i  and  N  is  the  number  of  combat  sorties.  Note  that  the 
stock  levels  computed  for  each  part  are  independent  of  the  stock  levels  computed 
for  other  parts. 

The  RBS-Part  model  is  analogous  to  the  readiness-based  spares  model  used 
by  ASO.  Given  a  stock  level  s,  the  probability  that  there  are  no  stockouts  over 
a  period  of  N  combat  sorties  is  taken  as  a  measure  of  the  readiness  provided  by 
the  stockpile.  Because  independence  is  assumed  in  this  model,  the  probability  of 
no  stockouts  may  be  written 


ns((i%, -■;,*)  . 

- _ 1 


An  inductive  procedure  is  used  to  build  a  stock  level  from  a  starting  position 
io  =  0.  At  each  step  in  the  process,  one  unit  of  the  item  offering  the  largest 
gain  in  “readiness”  per  unit  cost  is  added  to  the  stockpile.  More  precisely,  let  fn 
describe  the  stock  level  after  step  n.  Pick  M  to  satisfy 
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1 

(5)* 


n. 


II  *((£)i  4-  6jk,r„N)  -  n  B{(Z)JtrjtN) 
b= i  i=i  > 


^  (n  b{(  &n  )  j  4"  fijMil'jiN')  ni(w„'„jv) 


for  all  k  =  1,  •  •  •  ,np  where 


1  if  i  =  j 
0  otherwise 


Let  sn+i  =  sZ  4-  (£im>  •  •  • ,  6n,Ai)-  The  process  is  halted  when  a  specified  cost  goal 
is  reached. 

Both  of  the  above  spares  models  assume  independence  between  part  “fail¬ 
ures.”  This  assumption  does  not  seem  reasonable  when  applied  to  battle  damage. 
The  raw  output  of  the  NWC/NADC  modeling  effort  will  in  effect  be  stated  in 
terms  of  damage  incidents  and  must  be  condensed  into  part-level  damage  rates. 
Moreover,  it  seems  more  reasonable  to  assume  the  independence  of  damage  in¬ 
cidents  than  the  independence  of  part  damage.  This  motivates  the  following 
alternative  spares  models. 

The  DBS-Damage  Incident  model  is  a  direct  adaptation  of  the  DBS-Part 
model,  where  damage  incidents  rather  than  demands  for  parts  are  monitored.  In 
other  words,  each  damage  incident  corresponds  to  a  “bundle”  of  parts  needed  to 
repair  that  particular  damage.  Enough  bundles  of  parts  are  stocked  to  ensure 
that  demands  for  that  particular  bundle  are  met  with  at  least  a  probability  of 
90  percent.  More  precisely,  the  “stock”  level  of  damage-incident  i  is  the  first 
integer  k  where  B{k,(r‘p)i^N)  >  0.9.  The  damage-incident  “stock  levels”  Sd  are 
converted  to  part  stock  levels  i*  via  s  = 

In  an  analogous  manner,  the  RBS-Damage  Incident  model  is  adapted  from 
the  RBS-Part  model. 
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DAMAGE  RATE,  COST,  AND  PART-TO-DAMAGE- 
INCIDENT  DISTRIBUTIONS 


This  section  defines  the  damage-rate  distributions,  part-to-damage-incident 
distributions,  and  cost  distributions  used  in  the  analysis. 


F-14  Assemblies 

Initial  work  began  with  a  list  of  23  candidate  major  structural  assemblies  for 
ABDR  on  the  F-14  and  a  hypothetical  distribution  of  the  assemblies  among  30 
distinct  damage  incidents.  The  list  of  assemblies  and  costs,  listed  in  table  4,  were 
provided  by  NAVAIR.  The  part-to-damage-incident  table  is  given  in  table  5. 


Generalized  Part-to-Damage-Incident  Distributions 

Computer  limitations  would  not  allow  the  c<  nsideration  of  substantially  higher 
numbers  of  parts  and  damage  incidents;  hence,  23  parts  and  30  damage  incidents 
were  used  throughout  the  analysis.  However,  it  was  relatively  easy  to  consider 
more  general  part-to-damage-incident  distributions  than  the  hypothetical  F-14 
assemblies  distributions  shown  in  table  5. 

Given  a  fixed  number  r  €  (0, 1),  the  following  procedure  was  used  to  generate 
a  part-to-damage-incident  matrix  D.  For  each  fixed  damage  incident,  a  random 
number  was  drawn  for  each  part  j  to  determine  if  part  j  was  needed  to  repair 
the  damage  incident.  The  probability  that  part  j  was  needed  was  exactly  r.  If, 
at  the  end  of  this  process,  at  least  one  part  was  needed  to  repair  the  damage 
incident,  the  list  of  parts  just  selected  was  accepted.  If  no  parts  were  selected, 
the  (empty)  list  of  parts  was  rejected  and  the  selection  process  was  repeated. 

The  probability  distribution  function  corresponding  to  the  procedure  above 
is  given  by  the  following  function: 


t 

Prob{number  of  parts  included  >  x} 


—  F(r,x) 


1  -  B(x  —  l,r, 23) 
1  -5(0,  r,  23) 
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TABLE  4 


MAJOR  ASSEMBLIES  AND  COSTS 
FOR  CANDIDATE  F-14 


Cost 

Major  assembly _  (thousand  $) 


Outer  wing  panel  (left) 

617 

Outer  wing  panel  (right) 

617 

Aileron  (left) 

780 

Aileron  (right) 

780 

Leading  edge  flap  (left) 

28 

Leading  edge  flap  (right) 

28 

Trailing  edge  flap  (left) 

21 

Trailing  edge  flap  (right) 

21 

Horizontal  stabilizer  (left) 

445 

Horizontal  stabilizer  (right) 

445 

Vertical  tail  (left) 

106 

Vertical  tail  (right) 

132 

Rudder  (left) 

42 

Rudder  (right) 

42 

Main  landing  gear  (left) 

43 

Main  landing  gear  (right) 

109 

Engine  inlet  assembly  (left) 

30 

Engine  inlet  assembly  (right) 

30 

Canopy 

143 

Windscreen 

28 

Nose  landing  gear 

92 

Afterburner  assembly 

11 

Radome 

35 
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TABLE  5 


HYPOTHETICAL  F-14  DAMAGE  INCIDENTS 

Damage 

incident  number  Parts  requiring  replacement _ 

1-23  One  for  each  of  the  23  candidate  parts 

24  Outer  wing  panel  (left) 

Aileron  (left) 

Leading  edge  flap  (left) 

Trailing  edge  flap  (left) 

25  Outer  wing  panel  (right) 

Aileron  (right) 

Leading  edge  flap  (right) 

Trailing  edge  flap  (right) 

26  Canopy 
Windscreen 

27  Main  landing  gear  (left) 

Engine  inlet  assembly  (left) 

28  Main  landing  gear  (right) 

Engine  inlet  assembly  (right) 

29  Horizontal  stabilizer  (left) 

Vertical  tail  (left) 

Rudder  (left) 

Afterburner  assembly 

30  Horizontal  stabilizer  (right) 

Vertical  tail  (right) 

Rudder  (right) 

Afterburner  assembly 
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Each  of  the  23  candidate  parts  has  an  equal  probability  of  being  one  of  the 
x  parts  selected  for  a  damage  incident.  Three  different  values  of  the  parameter 
r  are  used  in  the  analysis.  The  parameter  value  r  =  0.023  was  selected  to 
match  with  the  hypothetical  F-14  assemblies  distribution  shown  in  table  5:  23 
percent  of  the  damage  incidents  include  two  or  more  parts.  The  parameter  value 
r  =  0.11  was  selected  so  the  probability  of  having  more  than  two  parts  included 
in  a  damage  incident  was  50  percent.  The  parameter  value  r  =  0.2  was  selected 
so  the  probability  of  having  more  than  four  parts  included  in  a  damage  incident 
was  50  percent.  On  average,  each  part  is  included  in  1.6  damage  incidents  when 
r  =  0.023,  3.5  damage  incidents  when  r  =  0.11,  and  6.0  damage  incidents  when 
r  =  0.2. 


Generalized  Cost  Distribution 


The  following  distribution  of  cost  was  used  to  generalize  the  cost  distribution: 


Prob 


=  < 


At 

0.8 

0.4i  +  0.6 


0  <  i  <  0.2 
0.2  <  t  <  0.5 
0.5  <  t  <  1 


This  distribution  was  chosen  to  match  reasonably  well  with  the  distributions 
of  major  aircraft  assembly  costs  provided  by  NAVAIR.  Figure  3  compares  the 
cost  distribution  function  with  the  actual  distribution  of  major  aircraft  assembly 
costs  for  the  A-6,  F-14,  and  F-18. 


Damage-Rate  Distributions 

A  two-parameter  distribution  of  predicted  and  actual  damage  rates  was  used 
in  this  analysis.  The  “maximum  DI  rate”  parameter  S  defined  the  range  of 
possible  damage  rates;  the  “confidence  factor”  C  defined  the  relationship  between 
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Prob  (CiC, 
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the  predicted  and  actual  damage  rates.  Precisely,  the  range  of  predicted  (rp)  and 
actual  (ra)  damage  rates  for  each  of  the  30  damage  incidents  actually  requiring 
a  part  was  uniformly  distributed  in  the  set 


{(rp,ra) :  0  <  rp  <  S  and  Crp  <ra  <  Crp  +  5(1  -  C )}  . 

A  parameter  value  C  =  0  corresponds  to  “no  confidence,”  a  value  C  =  0.5  corre¬ 
sponds  to  “some  confidence,”  and  a  value  C  =  1  corresponds  to  “full  confidence.” 

Figure  4  displays  pictorially  the  effect  of  the  confidence-factor  parameter  and 
the  maximum  DI  rate  parameter.  The  maximum  DI  rate  parameter  5  sets  the 
maximum  allowable  damage  rate.  For  example,  if  5  =  0.005,  the  probability  of 
experiencing  a  specific  damage  incident  on  a  sortie  is  no  more  than  0.005.  Note 
that  C  —  1  (full  confidence)  forces  the  predicted  and  actual  rates  to  be  equal,  and 
that  the  predicted  and  actual  rates  are  unrelated  when  C  =  0  (no  confidence). 
When  C  =  0.5  (some  confidence),  the  actual  damage  rate  and  predicted  damage 
rate  both  fall  within  an  interval  of  width  |. 
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Step  1 :  Pick  predicted  damage  rate  rp  in  interval  [0,  S]. 

Step  2:  Pick  actual  damage  rate  ra.  The  range  of  permitted  values  for  ra  is  shaded  below. 


No  confidence 
( C-0 ) 


Some  confidence 
(Cm  0.5) 


Full  confidence 
(C-1) 


o 


t 

r* 


S 


0 


t 


J 

S 


RG.  4:  EFFECT  OF  THE  CONFIDENCE-FACTOR  PARAMETER 
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APPENDIX  A 


DETAILED  RESULTS 


APPENDIX  A 


DETAILED  RESULTS 


This  appendix  contains  comparisons  of  the  four  ABDR  spares  requirements 
models  over  a  wide  range  of  parameter  and  distribution  values.  A  total  of  eight 
different  combinations  of  maximum  DI  rate,  part-to-DI  distribution,  and  cost 
distribution  were  considered  in  the  analysis,  as  listed  in  table  A-l.  The  term 
“series,”  defined  by  this  table,  is  used  to  refer  to  these  combinations.  Each  of  the 
tables  that  follow  compares  two  of  the  models  in  terms  of  the  statistic  “proba¬ 
bility  Model  A  performs  better  than  Model  B.”  More  precisely,  this  statistic  is 
Prob{CANNfl  -  CANNa  >  0}. 

The  relative  costs  of  the  stockpiles  associated  with  the  spares  models  must 
be  taken  into  account  when  comparing  the  spares  models.  The  cost  of  the  DBS- 
Damage  Incident  stockpile  is  always  higher  than  the  cost  of  the  corresponding 
DBS-Part  stockpile.  The  cost  associated  with  the  RBS  stockpile  is  indicated  at 
the  top  of  each  table.  For  most  of  the  comparisons,  the  cost  of  the  RBS  stockpile 
was  forced  to  remain  slightly  below  the  cost  goal:  the  RBS  algorithm  was  halted 
when  the  addition  of  the  next  part  would  exceed  the  cost  goal;  this  last  part  was 
not  added  to  the  stockpile.  This  rule  was  used  in  all  of  the  following  comparisons 
except  the  equal-cost  comparisons  of  the  RBS-Part  and  RBS-Damage  Incident 
models  (tables  A-8  and  A-9). 

In  the  equal-cost  comparisons  between  the  RBS-Part  and  RBS-Damage  In¬ 
cident  models,  the  cost  of  the  RBS  stockpiles  was  allowed  to  exceed  the  cost 
goal  slightly:  the  RBS  algorithm  was  halted  when  the  addition  of  the  next  part 
would  exceed  the  cost  goal;  this  last  part  was  added  to  the  stockpile.  This  case 
is  identified  in  the  following  tables  as:  “RBS  cost  is  slightly  higher  than  the  cost 
goal.”  The  “slightly  higher”  rule  causes  the  RBS-Damage  Incident  stockpile  cost 
to  be  slightly  higher  than  the  RBS-Part  cost  on  average.  The  other  rule  causes 
the  RBS-Part  stockpile  cost  to  exceed  the  RBS-Damage  Incident  cost  on  average. 
The  “slightly  higher”  rule  was  selected  for  the  equal-cost  comparisons  of  the  RBS 
models  because  it  gives  the  RBS-Damage  Incident  model  a  slight  advantage. 


A-l 


TABLE  A-l 


SERIES  COMPARED 


Series 

Max.  DI 

rate 

Part-to-DI 

dist. 

Cost 

dist. 

1 

0.0005 

F14  Assem. 

F14  Assem. 

2 

0.005 

F14  Assem. 

F14  Assem. 

3 

0.0005 

F14  Assem. 

Genl. 

4 

0.005 

F14  Assem. 

Genl. 

5 

0.0005 

Genl.,  r  =  0.023 

Genl. 

6 

0.005 

Genl.,  r  =  0.023 

Genl. 

7 

0.0005 

Genl.,  r  =  0.11 

Genl. 

8 

0.0005 

Genl.,  r  =  0.2 

Genl. 

9 

0.005 

Genl.,  t  =  0.2 

Genl. 

Tables  A-2  through  A-9  contain  comparisons  based  on  23  parts  distributed 
among  30  damage  incidents  and  a  support  period  of  3,000  sorties.  Discussions 
and  definitions  of  the  parameter  values  and  distributions  are  contained  in  the 
technical  section  of  this  paper.  The  accuracy  or  measurement  error  of  these 
statistics  is  also  discussed  in  the  technical  section  of  this  paper. 


TABLE  A-2 


DBS-P  VERSUS  RBS-P 
(Equal  cost) 


RBS  cost  goal  for  this  comparison 
RBS-P  cost:  DBS-P 


Conf. 

Series  factor 

1  No 
Some 
Full 

2  No 
Some 
Full 

3  No 
Some 
Full 

4  No 
Some 
Full 

5  No 
Some 
Full 

6  No 
Some 
Full 

7  No 
Some 
Full 

8  No 
Some 
Full 


Prob.  DBS-P 
is  better 

0.08 

0.09 

0.11 

0.08 

0.12 

0.20 

0.12 

0.13 

0.13 

0.12 

0.16 

0.26 

0.11 

0.11 

0.15 

0.13 

0.15 

0.24 

0.14 

0.17 

0.21 

0.14 

0.21 

0.19 


Prob.  RBS-P 
is  better 

0.53 

0.43 

0.38 

0.74 

0.72 

0.47 

0.44 

0.38 

0.33 

0.76 

0.68 

0.47 

0.46 

0.41 

0.32 

0.75 

0.70 

0.46 

0.49 

0.38 

0.32 

0.46 

0.36 

0.34 
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TABLE  A-3 


DBS-P  VERSUS  DBS-DI 
(DBS-P  cost  <  DBS-DI  cost) 

RBS  cost  goal  for  this  comparison 
Not  applicable 


Conf. 

Series  factor 

1  No 
Some 
Full 

2  No 
Some 
Full 

3  No 
Some 
Full 

4  No 
Some 
Full 

5  No 
Some 
Full 

6  No 
Some 
Full 

7  No 
Some 
Full 

8  No 
Some 
Full 


Prob.  DBS-P 
is  better 

0.02 

0.01 

0.01 

0.01 

0.00 

0.00 

0.02 

0.02 

0.01 

0.01 

0.00 

0.00 

0.01 

0.01 

0.00 

0.00 

0.00 

0.00 

0.01 

0.01 

0.00 

0.00 

0.00 

0.00 


Prob.  DBS-DI 
is  better 

0.34 

0.34 

0.40 

0.76 

0.73 

0.65 

0.30 

0.33 

0.36 

0.77 

0.76 

0.66 

0.33 

0.29 

0.32 

0.71 

0.61 

0.49 

0.58 

0.51 

0.51 

0.69 

0.62 

0.57 
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TABLE  A-4 


DBS-P  VERSUS  RBS-DI 
(DBS-P  cost  <  RBS-DI  cost) 

RBS  cost  goal  for  this  comparison 
RBS-DI  cost:  DBS-DI 


Conf.  Prob.  DBS-P 
Series  factor  is  better 


Prob.  RBS-DI 
is  better 


1 

No 

0.07 

0.60 

Some 

0.04 

0.52 

Full 

0.04 

0.47 

2 

No 

0.02 

0.89 

Some 

0.04 

0.86 

Full 

0.03 

0.65 

3 

No 

0.05 

0.55 

Some 

0.06 

0.51 

Full 

0.06 

0.47 

4 

No 

0.03 

0.88 

Some 

0.05 

0.83 

Full 

0.06 

0.65 

5 

No 

0.05 

0.53 

Some 

0.06 

0.48 

Full 

0.09 

0.40 

6 

No 

0.04 

0.87 

Some 

0.06 

0.80 

Full 

0.12 

0.55 

7 

No 

0.03 

0.66 

Some 

0.03 

0.54 

Full 

0.02 

0.51 

8 

No 

0.01 

0.71 

Some 

0.00 

0.62 

Full 

0.57 
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TABLE  A-5 


RBS-P  VERSUS  DBS-DI 
(RBS-P  cost  <  DBS-DI  cost) 

RBS  cost  goal  for  this  comparison 
RBS-P  cost:  DBS-P 


Conf.  Prob.  RBS-P 
Series  factor  is  better 


Prob.  DBS-DI 
is  better 


1 

No 

0.40 

0.18 

Some 

0.29 

0.22 

Full 

0.20 

0.23 

2 

No 

0.67 

0.25 

Some 

0.58 

0.25 

Full 

0.26 

0.38 

3 

No 

0.29 

0.20 

Some 

0.24 

0.25 

Full 

0.16 

0.24 

4 

No 

0.60 

0.23 

Some 

0.49 

0.30 

Full 

0.20 

0.40 

5 

No  ' 

0.29 

0.21 

Some 

0.28 

0.21 

Full 

0.17 

0.25 

6 

No 

0.49 

0.31 

Some 

0.46 

0.32 

Full 

0.23 

0.39 

7 

No 

0.18 

0.36 

Some 

0.11 

0.40 

Full 

0.06 

0.40 

8 

No 

0.06 

0.49 

Some 

0.03 

0.49 

Full 

0.02 

0.45 
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TABLE  A-8 


RBS-P  VERSUS  RBS-DI 
(RBS-P  cost  <  RBS-DI  cost) 

RBS  cost  goal  for  this  comparison 
RBS-P  cost:  DBS-P 
RBS-DI  cost:  DBS-DI 


Series 

Conf. 

factor 

Prob.  RBS-P 
is  better 

Prob.  RBS-DI 
is  better 

1 

No 

0.08 

0.27 

Some 

0.06 

0.27 

Full 

0.02 

0.26 

2 

No 

0.09 

0.66 

Some 

0.08 

0.61 

Full 

0.02 

0.51 

3 

No 

0.06 

0.31 

Some 

0.04 

0.33 

Full 

0.03 

0.29 

4 

No 

0.08 

0.67 

Some 

0.06 

0.66 

Full 

0.02 

0.52 

5 

No 

0.06 

0.29 

Some 

0.06 

0.24 

Full 

0.05 

0.19 

6 

No 

0.08 

0.62 

Some 

0.14 

0.50 

Full 

0.10 

0.37 

7 

No 

0.08 

0.45 

Some 

0.06 

0.42 

Full 

0.04 

0.40 

8 

No 

0.03 

0.53 

Some 

0.01 

0.51 

Full 

0.01 

0.45 

A-7 


TABLE  A-7 


DBS-DI  VERSUS  RBS-DI 
(Equal  cost) 


RBS  cost  goal  for  this  comparison 
RBS-DI  cost:  DBS-DI 


Conf.  Prob.  DBS-DI  Prob.  RBS-DI 
Series  factor  is  better  is  better 


1 

No 

0.14 

0.50 

Some 

0.13 

0.37 

Full 

0.13 

0.28 

2 

No 

0.13 

0.75 

Some 

0.12 

0.70 

Full 

0.17 

0.36 

3 

No 

0.11 

0.38 

Some 

0.11 

0.34 

Full 

0.13 

0.27 

4 

No 

0.12 

0.73 

Some 

0.15 

0.68 

Full 

0.18 

0.32 

5 

No 

0.11 

0.37 

Some 

0.12 

0.34 

Full 

0.19 

0.22 

6 

No 

0.12 

0.72 

Some 

0.16 

0.61 

Full 

0.24 

0.30 

7 

No 

0.13 

0.29 

Some 

0.11 

0.18 

Full 

0.10 

0.10 

8 

No 

0.10 

0.18 

Some 

0.09 

0.10 

Full 

0.05 

0.06 
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TABLE  A-8 


RBS-P  VERSUS  RBS-DI 
(Equal  cost,  “low”  cost  goal) 

RBS  cost  goal  for  this  comparison 
RBS-P  cost:  slightly  higher  than  DBS-P 
RBS-DI  cost:  slightly  higher  than  DBS-P 

Conf.  Prob.  RBS-P  Prob.  RBS-DI 
Series  factor  is  better  is  better 


1 

No 

0.21 

0.10 

Some 

0.17 

0.09 

Full 

0.13 

0.10 

2 

No 

0.33 

0.30 

Some 

0.28 

0.21 

Full 

0.22 

0.14 

3 

No 

0.17 

0.09 

Some 

0.17 

0.10 

Full 

0.13 

0.09 

4 

No 

0.30 

0.28 

Some 

0.26 

0.26 

Full 

0.23 

0.16 

5 

No 

0.20 

0.12 

Some 

0.25 

0.09 

Full 

0.24 

0.09 

6 

No 

0.31 

0.42 

Some 

0.38 

0.30 

Full 

0.40 

0.14 

7 

No 

0.41 

0.09 

Some 

0.32 

0.11 

Full 

0.31 

0.10 

8 

No 

0.39 

0.11 

Some 

0.38 

0.10 

Full 

0.31 

0.15 
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TABLE  A-9 


RBS-P  VERSUS  RBS-DI 
(Equal  cost,  “high”  cost  goal) 

RBS  cost  goal  for  this  comparison 
RBS-P  cost:  slightly  higher  than  DBS-DI 
RBS-DI  cost:  slightly  higher  than  DBS-DI 

Conf.  Prob.  RBS-P  Prob.  RBS-DI 
Series  factor  is  better  is  better 


1 

No 

0.15 

0.10 

Some 

0.14 

0.06 

Full 

0.08 

0.05 

2 

No 

0.33 

0.29 

Some 

0.29 

0.24 

Full 

0.12 

0.08 

3 

No 

0.18 

0.07 

Some 

0.14 

0.06 

Full 

0.10 

0.06 

4 

No 

0.27 

0.28 

Some 

0.27 

0.27 

Full 

0.15 

0.11 

5 

No 

0.25 

0.08 

Some 

0.23 

0.05 

Full 

0.18 

0.05 

6 

No 

0.33 

0.42 

Some 

0.46 

0.22 

Full 

0.32 

0.09 

7 

No 

0.31 

0.08 

Some 

0.25 

0.05 

Full 

0.14 

0.03 

8 

No 

0.23 

0.02 

Some 

0.15 

0.02 

Full 

0.07 

0.01 
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AN  EXAMPLE  COMPARING  RBS-PART 
AND  RB  S-D  AM  AGE  INCIDENT  MODELS 

Table  B-l  presents  a  more  detailed  example  of  the  reasons  for  the  excellent 
performance  of  the  RBS-Part  model  in  comparison  to  the  RBS-Damage  Incident 
model.  The  table  contains  estimates  made  under  the  “series  1”  conditions  of 
appendix  A  with  some  confidence  in  predicted  damage  rates.  In  particular,  the 
F-14  assemblies  cost  and  part-to-damage-incident  matrix  was  used.  These  pa¬ 
rameters,  taken  from  tables  4  and  5  of  the  main  text,  are  repeated  in  columns 
1  through  3  of  table  B-l.  The  other  columns  of  table  B-l  contain  estimates  of 
stockpile  characteristics  and  performance: 

•  Columns  4  through  6  of  table  B-l  contain  comparisons  of  the  stock  lev¬ 
els  produced  by  the  RBS-Damage  Incident  and  RBS-Part  models.  For 
example,  the  probability  that  part  1  is  stocked  to  a  lower  level  under  RBS- 
Damage  Incident  than  RBS-Part  was  estimated  to  be  0.28.  This  is  reported 
under  part  number  1  in  the  column  headed  “DI<P.” 

•  Columns  7  through  10  contain  the  probability  of  certain  types  of  stockout 
conditions.  For  example,  the  probability  of  a  stockout  of  part  1  in  the 
RBS-Part  stockpile  is  0.09.  This  is  reported  under  part  number  1  in  the 
column  labeled  “P.”  The  probability  of  a  stockout  for  part  1  in  the  RBS- 
Damage  Incident  stockpile  that  would  not  have  resulted  in  a  stockout  in 
the  RBS-Part  stockpile  is  0.04.  This  is  reported  under  part  1  in  the  column 
labeled  “DI,  not  P.” 

The  RBS  models  tend  to  trade  expensive  items  for  additional  protection  in 
inexpensive  items.  Intuitively,  the  RBS-Damage  Incident  model  will  find  damage 
incidents  26  through  28  attractive  because  of  the  low  cost  and  the  increase  in 
readiness  produced  by  adding  two  parts  to  the  stockpile  simultaneously.  Also 
intuitively,  the  RBS-Damage  Incident  model  will  not  find  parts  1  through  4  at¬ 
tractive  because  they  are  expensive  both  when  considered  separately  and  as  part 
of  damage  incidents  24  through  25.  By  this  intuitive  argument,  RBS-Damage 
Incident  will  tend  to  buy  fewer  of  parts  1  through  4  and  more  of  parts  15  through 
20  than  RBS-Part. 

The  intuitive  argument  above  agrees  with  the  estimates  presented  in  ta¬ 
ble  B-l.  Note  that  the  RBS-Damage  Incident  model  tends  to  stock  fewer  of 
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the  expensive  parts  and  parts  included  in  expensive  damage  incidents  (puts  1 
through  8)  than  the  RBS-Part  model,  and  to  stock  more  inexpensive  parts  and 
parts  included  in  inexpensive  damage  incidents  (particularly  parts  15  through 
20).  Consequently,  relative  to  the  RBS-Part  stockpile,  the  RBS-Damage  Inci¬ 
dent  stockpile  tends  to  experience  more  stockouts  for  expensive  items  and  fewer 
stockouts  for  the  inexpensive  items. 

The  RBS-Damage  Incident  stockpile  generally  does  better  than  the  RBS-Part 
stockpile  in  meeting  demands  for  inexpensive  items.  However,  the  RBS-Part 
stockpile  already  does  quite  well  in  meeting  demands  for  inexpensive  items.  In 
relative  terms,  there  is  little  difference  in  performance,  with  a  probability  of 
around  0.01  of  a  stockout  under  the  RBS-Part  model  that  would  be  covered 
under  the  RBS-Damage  Incident  model. 

Both  models  have  traded  off  expensive  parts  for  additional  coverage  of  in¬ 
expensive  items.  Their  performance  in  meeting  demands  for  expensive  parts  is 
relatively  poor.  However,  the  RBS-Part  model  performs  better  than  the  RBS- 
Damage  Incident  model  in  this  performance  category.  In  relative  terms,  the 
difference  in  performance  is  large,  with  up  to  a  0.04  probability  of  stockout  un¬ 
der  the  RBS-Damage  Incident  model  that  would  be  covered  by  the  RBS-Part 
model. 

In  summary,  the  RBS-Damage  Incident  model  appears  to  carry  the  tendency 
of  trading  off  expensive  items  for  inexpensive  items  a  bit  too  far.  Relative  to  the 
RBS-Part  model,  the  RBS-Damage  Incident  model  appears  to  experience  more 
stockouts  for  expensive  items.  While  it  is  true  that  the  RBS-Damage  Incident 
model  appears  to  experience  fewer  stockouts  for  inexpensive  items,  the  RBS- 
Part  model  also  does  well  in  meeting  demands  for  inexpensive  items.  The  overall 
effect  is  that  the  RBS-Part  model  generally  performs  slightly  better  than  the 
RBS-Damage  Incident  model. 


TABLE  B-l 


DETAILED  COMPARISON  OF  RBS-PART 
AND  RBS-DAMAGE  INCIDENT 


* 


» 


Comparison  of  stock  levels  Probability 

_ (probability  1 _  _ of  stockout 

Cost 


DI 

Part(s) 

(T$) 

DI<P 

O 

HH 

II 

T> 

DI>P 

DI 

P 

DI,  not  P 

P,  not  DI 

1 

1 

617 

0.28 

0.66 

0.06 

0.12 

0.09 

0.04 

0.00 

2 

2 

617 

0.27 

0.68 

0.05 

0.14 

0.11 

0.04 

0.01 

3 

3 

780 

0.20 

0.73 

0.07 

0.12 

0.11 

0.03 

0.01 

4 

4 

780 

0.21 

0.71 

0.08 

0.10 

0.08 

0.04 

0.01 

5 

5 

28 

0.62 

0.37 

0.01 

0.04 

0.01 

0.02 

0.00 

6 

6 

28 

0.61 

0.39 

0.01 

0.02 

0.01 

0.01 

0.00 

7 

7 

21 

0.66 

0.34 

0.01 

0.03 

0.01 

0.02 

0.00 

8 

8 

21 

0.68 

0.32 

0.00 

0.02 

0.01 

0.01 

0.00 

9 

9 

445 

0.02 

0.59 

0.39 

0.06 

0.07 

0.00 

0.02 

10 

10 

445 

0.04 

0.65 

0.31 

0.06 

0.08 

0.00 

0.02 

11 

11 

106 

0.16 

0.62 

0.22 

0.02 

0.01 

0.01 

0.00 

12 

12 

132 

0.12 

0.70 

0.19 

0.03 

0.03 

0.01 

0.00 

13 

13 

42 

0.25 

0.59 

0.15 

0.02 

0.01 

0.01 

0.00 

14 

14 

42 

0.23 

0.65 

0.11 

0.01 

0.01 

0.01 

0.00 

15 

15 

43 

0.00 

0.14 

0.86 

0.00 

0.01 

0.00 

0.01 

16 

16 

109 

0.00 

0.20 

0.80 

0.01 

0.02 

0.00 

0.01 

17 

17 

30 

0.01 

0.17 

0.82 

0.00 

0.01 

0.00 

0.01 

18 

18 

30 

0.02 

0.31 

0.67 

0.01 

0.00 

0.00 

0.00 

19 

19 

143 

0.01 

0.22 

0.77 

0.01 

0.02 

0.00 

0.01 

20 

20 

28 

0.02 

0.40 

0.58 

0.00 

0.00 

0.00 

0.00 

21 

21 

92 

0.09 

0.90 

0.01 

0.03 

0.03 

0.00 

0.00 

22 

22 

11 

0.37 

0.49 

0.14 

0.01 

0.00 

0.01 

0.00 

23 

23 

35 

0.08 

0.91 

0.01 

0.02 

0.01 

0.00 

0.00 

24  1,3, 5, 7  1,446 

25  2, 4, 6, 8  1,446 

26  19,20  171 

27  15,17  73 

28  16,18  139 

29  9,11,13,22  604 

30  10,12,14,22  630 
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main:  procedure  opt ions(main) ; 

/•«•  problem  definition  statements  •*•••*•••••••••••*••••••**/ 


declare 

printintermed 
print randno 
di dimens  ion 
partdimension 
sort  i.ed  intension 
samplesize 
scolefoctor 
conf factor 
seed 


character(80)  varying. 
character(80)  varying, 
fixed  binary, 
fixed  binary, 
fixed  binary, 
fixed  binary, 
float  binary, 
float  binary. 

fixed  binary,  /•  random  number  seed  •/ 


outf  i  le 
outf i le.name 
1  outrec, 

2  out_i 
2  out.cannl 
2  out_cann2 
2  out_cann3 
2  out_cann4 
2  out.costl 
2  out_cost2 
2  out_cost3 
2  out_cost4 


file, 

char(80)  varying, 

picture  ‘zz9‘, 
picture  'zzzz9*. 
picture  'zzzz9', 
picture  ‘zzzz9', 
picture  'zzzz9', 
picture  'zzzzzzz9v.9‘ . 
picture  'zzzzzzz9v.9‘ , 
picture  *zzzzzzz9v.9’ , 
picture  ‘ zzzzzzz9v.9’ . 


p_accept_part 

cost_phigh 

cost_low_l 

cost_low_u 

cost_high_l 

cost_high_u 


float  binary, 

float  binary, 
float  binary, 
float  binary, 
float  binary, 
float  binary; 


I 


I 


get  I ist (printintermed) 

get  skip  I ist(printrandno) 
get  skip  I ist(didimension) 
get  skip  I ist(partdimenslon) 
get  skip  I i st (sort led imens ion) 
get  skip  I ist(samplesize) 
get  skip  I ist(scalefactor) 
get  skip  I ist(conf factor) 
get  skip  I ist (seed) 
get  skip  I ist (outf i le_nome) 
get  skip  I lstfp_aecept_port) 
get  skip  I ist(cost_phigh) 
get  skip  I lst(cost_low_l ) 
get  skip  I ist(cost_low_u) 
get  skip  I ist(cost_hfgh_l) 
get  skip  I ist(cost_high_u) 


opt i ons( prompt ( ' pr i nt i ntermed  (Y/N) ’ ) ) ; 

opt  ions (prompt ( 'print randno  (Y/N) ’ )) ; 

opt  ions (prompt ( 'di dimens  ion’ )) ; 

opt ions(prompt( 'partdimension' )) ; 

opt  ions (prompt ( 'sort  led imens ion' )) ; 

opt i ons( p rompt ( ' samp  I es i ze ' ) ) ; 

opt ions(prompt( ’ sealefactor ' ) ) ; 

opt ions(prompt( 'conf factor ’)) ; 

opt ions(prompt( 'seed  (-1  for  svtem  seed)')); 

opt ions(prompt( 'outf i le.name’ )) ; 

opt ions(prompt( 'p_occept_port ■ )) ; 

opt  ions (prompt ( ' cost_phi gh ' ) ) ; 

opt i ons(prompt ( ' cost_l ow_l ’)); 

opt ions(prompt( 'cost.low.u')) ; 

opt ions(prompt( 'cost_high_l ’)); 

opt ions(prompt( 'cost.high.u' )) ; 


put  page  I ist ('echo  problem  definition  variables’); 
put  skip  I ist (did imens ion, 'did Imens ion’ ) ; 
put  skip  I ist(partdimension, 'partdimension'); 
put  skip  I ist(sort iedimension, 'sort iedimension’ ) ; 
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put  skip  I ist(samplesi ze, ’ samples! ze' ) ;  t 

put  skip  I ist (seal sf actor . 'scalef actor ' ) ; 

put  skip  I ist (conf factor, 'conf factor*): 

put  skip  I i st (out f i I e.name, 'out f i le_name ' ) ; 

put  skip  I ist(p_occept_port , 'p_accopt_part ‘ ) ; 

put  skip  I ist(cost_phigh, 'cost.phigh' ) ;  f 

put  skip  I ist (cost_low_l , 'cost_l ow_l ’ ) ; 

put  skip  I i st (cos t_l ow_u , ‘ cost_l ow_u ’ ) ; 

put  skip  I ist(cost_high_l , 'cost_high_l '); 

put  skip  I i st (cost_hi gh_u , ’cost_high_u' ) : 

beg i n ; 

open  f i I e(out f i I e)  t i t le(outf i te_nane)  output  record 

envi ronment(f i xed_l ength_records ,max imum_record_si ze(63)) ; 


/•••  declaration  of  globals  variob 


'/ 


declare 

di (part  d  imens  i  on ,  d  i  d  inn 
c(partdimension) 

ra(didimension) 
rp(di dimension) 
rk(sort iedimension) 

rnset  external  entry 

rnget  external  entry 

rnunf  external  entry 


ion)  fixed  binary.  /• 
float  binary,  /• 

float  binary(113), 
float  binary(l13), 
float  binary(l13), 

(fixed  binary).  /• 
(fixed  binary),  /• 
returns(f (oat  binary);  /• 


part-di  cross  reference  •/ 
cost  vector  •/ 

/•  actual  rates  •/ 

/•  predicted  rates  •/ 

/•  damage  outcome  •/ 

set  seed  •/ 
get  seed  •/ 
random  number  •/ 


if  seed  >  0  then  call  rnset(seed); 
call  rnget(seed): 
put  page; 

put  skip  I ist(' random  number  seed:  '.seed); 

/•••••••  main  routine  . . . . •••••*••••***/ 

/••••«•  monte  carlo  integration  •••••••••••••••••••••••••••••••••/ 

declare 


i 

fixed  binary. 

s1(partdimension] 

1  fixed  binary. 

/• 

spares  level  model  1  •/ 

s2(part dimens  ion] 

1  fixed  binary. 

/• 

spares  level  model  2  •/ 

s3(partd (mens  ion] 

1  fixed  binary. 

/• 

spares  level  model  3  •/ 

s4(partdimension] 

1  fixed  binary. 

/• 

spares  level  model  4  •/ 

k(didimension) 

fixed  binary. 

/• 

damage  outcome  •/ 

cannl 

fixed  binary. 

/• 

conn  ac  for  model  1  •/ 

cann2 

fixed  binary. 

/• 

conn  ac  for  model  2  •/ 

cann3 

fixed  binary. 

/• 

conn  ac  for  model  3  */ 

cann4 

fixed  binary. 

/• 

conn  ac  for  model  4  •/ 

costl 

float  binary. 

/• 

cost  of  model  1  parts  •/ 

cost2 

float  binary. 

/• 

cost  of  model  2  parts  •/ 

cost3 

float  binary. 

/• 

cost  of  model  3  parts  •/ 

cost4 

float  binary; 

/• 

cost  of  model  4  parts  •/ 

print intermed  «  'V  then  put  skip(3)  1 i st ( * 

intermediate  results’); 
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do  i  •  1  to  aompleaize; 


call  getcoatvect; 

call  getdimatrix; 

call  getdamagavect ; 

call  a_aao(a1); 

call  computecoat(a1 .cost  1 ) ; 

call  a_rbs(a2^cost1) ; 

call  computeco3t(a2,coat2) ; 

call  a_aaodi(a3); 

call  computecoat(a3,coat3) ; 

call  s_rbadi (a4,coat3) : 

call  conputecoat(a4,co3t4); 

coll  computak(k) ; 

call  computecann(k,a1 .cannl); 

call  conputecann(k,a2,cann2); 

call  computecann(k . a3,cann3) ; 

call  computecann(k,a4,cann4) ; 


if  print interned  -  *Y* 
than 
do; 

put  skip(2) ; 
put  skip  I iat( * i :  ' .  i); 
put  akip  I i at ( * al :  *,a1); 
put  akip  Iiat(‘a2:  *.a2); 
put  akip  Iiat('a3:  ' ,a3); 
put  akip  Iiat(’a4:  • ,a4); 
put  akip  liat(’coat1:  ’.coatl); 
put  akip  liat('eoat2:  ’,coat2); 
put  akip  liat('coat3:  ' ,coat3); 
put  akip  liat(’coat4:  ’ ,coat4); 
put  akip  I iat(*k:  *,k); 
put  akip  I i at( ’cannl : ’ , cannl ) ; 
put  akip  I iat('conn2: * ,cann2); 
put  akip  I iat('conn3: ' ,eann3); 
put  akip  I iat(’cann4: • ,cann4); 

and; 

out_i  -  i ; 

out.cannl  -  cannl ; 

out_cann2  “  cann2; 

out_cann3  -  cann3; 

out_cann4  ■  cann4; 

out.coatl  -  coatl; 

out_coat2  -  coat 2; 

out_coat3  “  coat3; 

out_coat4  ”  eoat4; 

write  f i le(outf i la)  f ron(outrac) ; 

and; 

stop; 

/ . * . / 

computak;  procadura(k) ; 
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declare 

k(.) 


fixed  binary,  /«  damage  outcome  •/ 

fixed  binary, 
fixed  binary, 
float  binary; 


do  i  >  1  to  didimension; 
k(i)  -  0; 

end; 

do  i  ■  t  to  sort iedimens ion; 
r  -  rk( i ) ; 

j  -  0; 

do  whi le(r  >  0); 

j  -  j  +  1; 

r  -  r  -  ra( j ) ; 

end; 

k(j)  -  k(j)  +  1; 

end; 

end  computek; 

/••«••  compute  number  of  cannibalized  aircraft 

computecann:  procedure(k.s.nc) ; 

declare 


3(0 

fixed  binary. 

k(0 

fixed  binary, 

nc 

fixed  binory. 

i 

fixed  binary. 

i 

fixed  binary. 

dp 

fixed  binary; 

/•  spares  level  •/ 

/•  damage  outcome  •/ 

/•  number  of  conned  oc  •/ 


nc  ■  0; 

do  i  «  1  to  partdimension; 
dp  -  0; 

do  j  «  1  to  didimension; 
dp  - dp  +  k(j)edi(i , j); 

end; 

nc  “  max(nc,dp-s( i )) ; 

end; 

end  computecann; 


getdimatrix:  procedure; 
declare 

acceptdi  bit(1), 

rvect(partdimeneion,2;didimension)  float  binary, 

i  fixed  binary, 

j  fixed  binary; 


/•  damage-part  cross  reference  di (port .damage  incident)  •/ 
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/•  logic  assumes  damage  incident  1  is  the  ’no  damage'  event  •/ 


do  i  -  1  to  portdimension; 
di(i.l)  -  0; 

end; 

do  j  ■  2  to  didimension; 
acceptdi  -  '0’b; 
do  whi le(tacceptdi ) ; 
do  i  •  t  to  portdimension; 
rvect(i ,  j)  *■  rnunf(); 
if  rvect(i.j)  <  p_accept_part 
then 
do; 

di  (  i .  j )  ■  1 : 

acceptdi  ”  *  1 ’ b ; 
end; 

else  di(i,j)  -  0; 

end; 

end; 

end; 

if  printrandno  ■  *Y‘ 
then 

do  i  -  1  to  partdimension; 

j  “  1: 

put  skip  list(>i.j.*.di(i,j)’,i,j,di(i,j)); 
do  j  -  2  to  didimension; 

put  skip  I i st ( ’ i . j . rveet( i ,j).di(i.j)',i.j. rvect(i , j ) ,di ( i , j )) ; 

end; 

end; 

end  getdimatrix; 

/ . ••*•/ 

getcostvect:  procedure; 

/•  cost  vector  c(port)  •/ 
declare 

rvect(2. portdimension)  float  binary, 

I  fixed  binary; 

do  i  -  1  to  portdimension; 
rvect(l.i)  -  rnunfQ; 
rvect(2,i)  -  rnunfQ: 

end; 

do  I  *  1  to  partdimension; 

if  rvect(t.i)  <  cost.phigh 
then 

c(i)  ■  cost_high_l  +  (cost_high_u-cost_high_l )*rvect(2, i) ; 

else 

c(i)  ■  cost_low_l  +  (eost_low_u-cost_low_l )*rvect(2, i ) ; 

end; 

I f  printrandno  ■  ’Y* 
then 
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do  i  -  1  to  partdimension; 

put  skip  I ist( ' i . r(*. i ) ,cost( i ) :  * . i ,rvect(1 . i),rvect(2, i).c(i)) 

snd; 

«nd  getcostvect; 

/••*•••  random  damage  rats  vector  generator  •••••••••••/ 

ge t damage vec t :  procedure; 

/•  logic  assumes  damage  incident  1  corresponds  to  no  damage  •/ 
declare 

rvect(didimension-1 )  float  binary, 

i  fixed  binary; 

/•  rp  •/ 

do  i  ■  1  to  di dimens  ion-1 ; 
rvect(i)  ■  rnunf(); 

end; 

rp(1 )  -  Is 

do  i  ■  2  to  di dimens  ion; 

rp(i)  ”  sea  I ef octor* rvec t ( i-1 ) ; 
rp(1)  -  rp(1)  -  rp( i ) ; 

end; 

If  printrandno  “  *Y* 
then 

do; 

put  skip  I ist ( ’ rp(1 ) ’ , rp(1 )) ; 
do  i  >  2  to  di dimension; 

put  skip  I ist( * i , r, rp:  * , i . rvect( i-1) , rp( i )) ; 

end; 

end; 


/•  ra  “  (1-conf factor)*ra  +  conf factor*rp)  •/ 

do  !  ■  I  to  di dimens  ion-1 ; 
rvect(i)  ■  rnunf(); 

end; 

ra(1)  -  1; 

do  i  ■  2  to  di dimens  ion; 

ra(i)  •  scalefactor*rvect( i-1 ) ; 

ro(i)  “  (1-conf factor)*ra( i )  +  conf f actor*rp( I ) ; 

ra(1)  -  ro(1 )  -  ro(i); 

end; 

if  printrandno  ■  *Y* 
then 

do; 

put  skip  I ist( ’ ra(1 ) ’ , ra(1 )) ; 

do  I  -  2  to  didimenslon; 

put  skip  I i st ( *  I .r.ro:  • , i , rvect( i-1 ) , ro( i )) ; 

end; 

end; 
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/•  «•*  •/ 


do  i  *  1  to  sort iedimension; 
rk(i)  -  rnunf(); 

and; 

if  printrondno  ■  *Y* 
then 

do  i  a  1  to  sort iedimension; 

put  skip  list('i.rk:  * » i . rk( » ) ) ; 

and; 

and  getdomagevect ; 


/•  oso  spares  policy  implementat ion  •/ 
s_oso:  procedure(s) ; 

Xreplace  safetyfactor  by  0.9; 


declare 

»(•) 


fixed  binary,  /•  spares  level  */ 

float  binary(113),  /•  part  damage  rates  •/ 


r_par  t  (pa  rtd  intension) 

pn  fixed  binary, 
j  fixed  binary, 
nm  float  bi nary( 113), 
bb  float  b inary(llJ) . 
tt  float  binary(113); 

/•  build  up  damage  rates  for  parts  •/ 

do  pn  m  1  to  partdimension; 
r_part(pn)  -  0; 
do  j  -  1  to  didimension; 

r_part(pn)  -  r_part(pn)  +  di (pn , j )• rp( j ) ; 

end; 

end; 

/•  built  aso  stock  level  •/ 

do  pn  -  1  to  partdimension; 
f_port(pn)  t-  1 
then 

do; 

mm  -  r_part(pn)/(1-r_part(pn)); 
bb  “  (1-r_part(pn))»«sortied(mension; 
tt  -  bb; 
s(pn)  -  0; 

do  while(bb  <  safetyfactor); 
s(pn)  -  s(pn)  +  1; 

tt  ”  tt*mm*(sort  iedimenaion-t1-s(pn))/s(pn)  ; 
bb  -  bb  +  tt; 

end; 

end; 

else 

do; 
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put  skip  I iat ('trivial  cos*. . . * . r_part) ; 
stop; 


•nd; 


tnd; 

•nd  s_aso;  f 

. . . . . 


/•  aso  damage  incident  spores  policy  implementation  •/ 
s.asodi :  procedure(s) ; 

Xreplace  safetyfactor  by  0.9; 
declare 

•(•)  fixed  binary,  /»  spores  level  •/ 


dis(2:di dimens  ion) 

fixed 

binary, 

din 

fixed 

binary. 

j 

fixed 

binary. 

mm 

float 

binary(113). 

bb 

float 

binary! 113) , 

tt 

float 

binary(113) ; 

/•  built  aso  damage  incident  stock  level  •/ 

/•  assumes  damage  incident  1  is  the  'no  damage’  event  •/ 

do  din  >  2  to  didimension; 
if  rp(din)  t«  1 
then 

do; 

mm  "  rp(din)/(1-rp(din)); 

bb  «  (1-rp(din))»*sort iedimension; 

tt  •  bb; 

dis(din)  -  0; 

do  while(bb  <  safetyfactor); 
dis(din)  ■  dis(din)  +  1; 

tt  “  tt»mm*(sort iedimension+1-di s(din))/di s(di n) ; 
bb  -  bb  +  tt; 

end; 

end; 

else 

do; 

put  skip  I ist('trivial  case — *.rp); 
stop; 

•nd; 

•nd; 

do  j  ■  1  to  partdimension; 

do  din  >  2  to  didimension; 

»(j)  ■  »(j)  +  di ( j ,din)*dis(din) ; 

•nd; 

•nd; 
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•nd  s.osodi ; 


/•  rba  aparaa  policy  implementation  •/ 
a_rba:  procedures, costgoal ) ; 


declare 

»(•) 

cost goo  I 


r_port (part dimension) 


fixed  binary, 
float  binary, 

float  binary(113). 


j  fixed  binary, 

im(partd intension)  float  binary(113). 

bb(partdimension)  float  binary(113), 

tt(partdimension)  float  binary(113), 

pn  fixed  binary. 

bestpn  fixed  binary, 

bestratio  float  binary(113), 

currdelta  float  binory(113), 

currratio  float  binary(113), 

cost  float  binary; 

/•  initial ise  •/ 

cost  »  0; 

/•  build  up  damage  rates  for  parts  •/ 

do  pn  »  1  to  partdimension; 
r_port(pn)  »  0; 
do  j  -  1  to  didimension; 

r_part(pn)  -  r_part(pn)  +  di(pn. j)erp(j); 

end; 

end; 

do  pn  m  t  to  partdimension; 
s(pn)  -  0; 

If  r_part(pn)  1 
then 

do; 

mm(pn)  -  r_port(pn)/(1-r_port(pn)) ; 
bb(pn)  -  (1-r_port(pn))»»sort iedime 
tt(pn)  -  bb(pn) ; 


mens i on ; 


put  skip  I i st ( ’Trivial  problem:  ’,r_part) 
stop; 


/•  build  rbs  level  •/ 
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do  whil«(eost  <  costgool); 

currdolto  -  t t (1 )«mm(1 )*(sort iedimensi on-s( 1 ))/(o(1 )+1 ) ; 

b«9' rin  «  1 ; 

bast  ratio  -  currdol to/(c(1 )«bb(1 )) : 
do  pn  -  2  to  partdimension; 

currdolto  *  t t (pn)*mm(pn)»(oort ied imens i on-s(pn) )/(o(pn)+1 ) ; 
currrotio  «  currdol ta/(c(pn)«bb(pn)) ; 
if  currrotio  >  bestratio 
thon 

do; 

beotpn  ■  pn; 
bootratio  »  currrotio; 

ond; 


ond; 

coot  -  coot  +  c(bostpn); 
if  coot  <•  costgool 
thon 
do; 

tt(bootpn)  *  t t (boo tpn)*mm(baotpn)*( sort i ed i mens i on-s( boo tpn))/(o( boo tpn)+1) ; 
bb(bootpn)  ■  bb(bestpn)  +  tt(bostpn); 
o(bootpn)  «  o(bootpn)  +  1; 

ond; 

ond; 


id  o_rbo; 


T 


r 


/ . . . . 

. . 

/•  rbo  damage  incident  policy  implomontat ion  •/ 

o.rbodi:  procoduro(o,cootgoal ) ; 


declare 

o(*)  fixed  binary, 

cootgoal  float  binary. 


dio(2:didimenoion) 


mm(2:di dimension) 
bb(2 : d i d imons i on) 
tt(2:didimension) 
din 

bostdin 

bootratio 

currdolto 

currrotio 

coot 

edi (2:di dimono  ion) 
/•  initialize  •/ 
coot  -  0; 

/•  initial izo  edi  •/ 


fixed  binary, 

fixed  binary, 
float  binary(113), 
float  binary(113), 
float  binary(H3)i 
fixed  binary, 

fixed  binary, 
float  binary(113), 
float  binary(113), 
float  b!nary(113), 
float  binary, 
float  binary; 


do  din  -  2  to  didimonoion; 


/ 


i 


\ 
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T 


1 


1 


T 


cdi(din)  -  0; 

do  j  ■  1  to  portdimenaion; 

cdi(din)  -  cdi(din)  +  di ( j , d i n)*c( j ) ; 

•nd; 

•nd; 


/•  aaaumea  damage  incident  1  ia  the  'no  damage’  event  •/ 

do  din  ■  2  to  didimenaion; 
dia(din)  -  0; 
if  rp(din)  t«  1 
then 

do; 

mm(din)  -  rp(din)/(1-rp(dln)) ; 
bb(din)  »  (1-rp(din))**aort iedimenaion; 
tt(din)  ■  bb(din); 

end; 

elae 

do; 

put  akip  I iat(’Trivlol  problem:  a.rp); 
atop; 

end; 

end; 


/•  build  rbadi  level  •/ 


do  while(coat  <  coatgool); 

currdel to  -  t t(2)*mm(2)*(aort iedimenaion-di s(2))/(di a(2)+1 ) ; 
beatdin  -  2; 

beatrotio  -  currdel to/(cdi (2)*bb(2)) ; 
do  din  -  2  to  didimenaion; 

currdel to  -  tt(din)*mm(din)«(aort iedimenaion-di a(din))/(dia(din)+1); 
currrotio  -  currde I to/(cdi (di n)*bb(di n)) ; 
if  currrotio  >  beatrotio 
then 


do; 

beatdin  »  din; 
beatrotio  -  currrotio; 

end; 

end; 

coat  ■  coat  +  cdi(beatdin); 
if  coat  <-  coatgool 
then 
do; 

tt (beatdin)  -  tt(beatdin)*mm(beatdin)*(aort iedimenaion-di a(beatd in) )/(dis(beatdi n)+1 ) ; 
bb(beatdin)  •  bb(beatdin)  +  tt(beatdin); 
dia(beatdin)  ”  dia(beatdin)  +  1; 

end; 

end; 


do  j  ■  1  to  portdimenaion; 

•(j)  -  «; 

do  din  -  2  to  didimenaion; 

•(j)  -  a(j)  +  di(j.din)edia(din); 

end; 

end; 


end  a.rbadi ; 
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t 


. . 

computecost:  procedural, coat)  ; 
declare 

i  fixed  binary, 

e(«)  fixed  binary, 

cost  float  binary; 

/•  compute  cost  */ 
cost  “  0; 

do  i  ■  I  to  portdimenslon; 
cost  ■  cost  -f  s(i)*c(i); 

end; 

end  computecost; 

/ . . . . 

end;  /•  begin  •/ 
end  main; 


I 


1 
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